CONSTRUCTING AN EXPANDING METRIC 
FOR DYNAMICAL SYSTEMS IN ONE COMPLEX VARIABLE 

o 

O ' SUZANNE LYNCH HRUSKA 

Abstract. We describe a rigorous computer algorithm for attempting to con- 
struct an explicit, discretized metric for which a polynomial map /: C — ♦ C 
is expansive on a neighborhood of the Julia set, J . We show construction of 
such a metric proves the map is hyperbolic. We also examine the question of 
whether the algorithm can be improved, and the related question of how to 
build a metric close to euclidean. Finally, we give several examples generated 
r^ • with our implementation of this algorithm. 
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1. Introduction 

Our main concern in this paper is to develop and use rigorous computer in- 
vestigations to study the dynamics of polynomial maps of C of degree d > 1; for 
^ ■ example, f{z) = z^ -\-c. For complex polynomial maps, the invariant set of interest 

(y^ \ is the Julia set, J . The Julia set can be defined as the topological boundary of the 

^D ■ set of points with bounded orbits. Intuitively, J is precisely where the chaotic dy- 

^^ \ namics occurs. For example, / is topologically transitive on J; also, J is non-empty 

and perfect (see for instance |TT]'). 
—J. \ The hyperbolic polynomials are a large class of maps with chaotic dynamics, but 

f^ ■ whose stability properties make them amenable to computer study. A polynomial 

map / is called hyperbolic, or expansive, if / is uniformly expanding on some neigh- 
borhood of J, with respect to some riemannian metric. Uniform expansion forces 
some dynamical rigidity. For example, hyperbolicity of a polynomial map / implies 
structural stability, i.e., in a neighborhood of / in parameter space, the dynamical 
behavior is of constant topological conjugacy type. Further, if / is hyperbolic then 
the orbit of every point in the complement of J tends to either some attracting 
periodic orbit or infinity. In fact, a polynomial map / is hyperbolic if and only if 
jrt ' the orbit of every critical point of / tends to an attracting periodic orbit or infinity. 

Thus the fate of the critical points provides a straightforward test for hyperbolicity. 
In this paper we develop an alternate test for hyperbolicity which produces more 
explicit information about the dynamics of any given map. 

We begin with the work of [^ as a foundation. There we described a rigorous 
algorithm (and its implementation) for constructing a neighborhood 13 of J , and 
a graph F which models the dynamics of / on 23. In this paper, we further our 
program by developing a rigorous algorithm for attempting to construct a metric in 
which a given / is expansive, by some uniform factor L > 1, on the neighborhood !B. 
We show that successful construction of such a metric proves that / is hyperbolic. 
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2 S. L. HRUSKA 

In addition, we analyze both some limitations of and possible improvements to 
our basic algorithm. We show, via counterexamples, that a significantly simpler 
algorithm will not suffice. On the other hand, we do present an enhancement to 
our algorithm, to build a metric closer to euclidean. 

One motivation for improving the algorithm and the metric which it builds is 
that in j^, we use this one dimensional algorithm as part of a computer-assisted 
proof of hyperbolicity of polynomial diffeomorphisms of C^, Hcnon mappings in 
particular. The one dimensional algorithm is used in such a way that that the 
closer to euclidean we can build the metric, the more likely we are to succeed at 
proving hyperbolicity in C^. Thus improving the current methods could lead to 
more interesting examples of proven hyperbolic polynomial diffeomorphisms of C^. 

We implemented the algorithms of this paper in a program Hypatia, ^ and at 
the conclusion of the paper, we give several examples of metrics constructed which 
establish hyperbolicity of some polynomial maps of degrees two and three. 

1.1. Statement of main results. In |^ we described the 60a: chain construction, 
for building a directed graph F representing / on some neighborhood !B of J. A 
similar approach is in the body of work described in the survey jE] (cf jSUHl^l^l)' 
In this paper we need only: 

Definition 1.1. Let F = (V, £) be a directed graph, with vertex set V — V{T) ~ 
{^k}k=iy ^ finite collection of closed boxes in C, having disjoint interiors, and such 
that the union of the boxes S = 'B(F) = U^^^Bk contains J. Suppose there is a 
(5 > such that F contains an edge from Bk to Bj if the image f{Bk) intersects a 
(5-neighborhood of i?j , i.e., 

eD{{k,j):f{Bk)ny^{B,,6)^^}. 

Further, assume F is strongly connected, i.e., for each pair of vertices Bi^, Bj, there 
is a path in F from Bk to Bj, and vice- versa. Then we call F a 60a; chain model of 
f on J and 3 a box .Julia set. 

For the remainder of this paper, F will denote a box chain model of / on J, for 
some polynomial map / : C ^ C of degree greater than one. 

Remark 1 . A box chain model of / on J satisfies the definition of a symbolic image 
of f, given by Osipenko in |15j . 

Given a box chain model F. our first goal in this paper is to define a version of 
hyperbolicity for F which is checkable by computer, and which implies hyperbolicity 
for /. Our discrete version of hyperbolicity is parallel to the following standard 
definition, from [111 : a polynomial / : C — » C is called hyperbolic if there exists 
a riemannian metric fj., on some neighborhood "B of J, and an expansion constant 
A > 1, such that the derivative D^f at every point z in Ti satisfies \Dzf['v)\ > 
A |v| . for every vector v in the tangent space T^C 

First, we introduce the following discrete version of a riemannian metric. 

Definition 1.2. Let V = {Bk}k=i be a finite collection of closed boxes in C, with 
disjoint interiors. Let $: (LJ^j^i?^) -^ IR+ be a box constant function, i.e., for some 
set of positive constants {fk}^^i, we have $(2;) = ipk H z G B^. Then <& times 
euclidean induces a metric on TV = U^^-fTBk. That is. for any Bk G V, z G Bk, 

To obtain a copy of this unix program, write to the author. 
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and V e TzBk, let |v|j. := ipk |v|. Let |-|^ be the induced norm on U^^^TBk- We 
call $ a box metric^ and say ipk is a handicap for _Bfe. 

Note that if 2 G Sfe n Bj, then |v|^ = |v|^ if we are considering v e T^Bk, but 
l"^!* ~ 1^1 7" '^^ ^^ ^^^ considering v G T^Bj. We justify the choice of the name 
"handicap" a few paragraphs below. 

Now wc present our discretized version of hyperbolicity. 

Definition 1.3. Call T = (V, £) box expansive if there exists a box metric $ on 

V and a 602; expansion constant L > 1 such that for all (fc,j) e £, z e Sfc, and 

V G T.Bfc, we have |i?./(v)|^. > L |v|^. 

The following theorem is a first piece of evidence that this definition is useful. 

Theorem 1.4. Suppose there exists a box metric $ and an L > 1 for which T is 
box expansive. Then f is hyperbolic. 

In particular, there exists a smooth function p: 23 (F) —t R+, which defines a 
riemannian metric on TTi (by p times euclidean) , such that for all z g 23, and 

V e r,S, we have \D,f{w)\p > L |v|^ . 

The key step in the proof of this theorem is to smooth out the box metric $ in a 
small neighborhood of the box boundaries, using a partition of unity argument and 
the "edge overlap" factor 6 from Definition ll.il The proof is given in Section |21 

We develop an algorithm for attempting to build a box metric for which a given 
r is box expansive, called the Handicap Hedging Algorithm, in Section 2J Then we 
describe how any outcome of the algorithm gives useful dynamical information, in 
Sectional For example, we obtain: 

Theorem 1.5. Given F and L > 0, the Handicap Hedging Algorithm either 

(1) constructs a box metric for which F is box expansive by L, or 

(2) produces a cycle of F which is an obstruction, showing there exists no box 
metric for which F is box expansive by L. 

Note that showing box expansion by L = 1 may be instructive, but it is not 
enough to prove hyperbolicity. Thus, if any algorithm, for example, the Handicap 
Hedging Algorithm, builds a box metric for which some F is box expansive by some 
L > 1, then / is proven hyperbolic. 

1.2. Secondary results and discussion of the approach. A natural question 
one might ask is whether there is a simple algorithm for showing box expansion, 
without constructing a box metric; for example, by examine the multipliers along 
the simple cycles in the graph. (A cycle in a graph is called simple if is it composed 
of distinct vertices.) In fact, examining cycle multipliers is the key idea to a better 
understanding of box expansion and box metrics. We would like to thank Clark 
Robinson for asking this question. 

Definition 1.6. Call Afe = min{|/'(z)| : z E Bk} the multiplier of Bk. If Bq -^ 
. . . ^ -B„_i -^ Bn = Bq is an n-cyclc of boxes in F, then the cycle multiplier is 
(Ao • • • A„_i), and the average cycle multiplier is (Aq • • • A„_i)^'". 

Along the way to proving Theorem 11.51 we establish the following characteri- 
zations of box expansion, which are independent of any algorithm used to find an 
expanded box metric. 
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Proposition 1.7. Let L be the minimum average cycle multiplier over all simple 
cycles in the graph T . Then for any L > 0, there exists a box metric for which T is 
box expansive by L if and only if L < L. 

Corollary 1.8. Let M be the minimum cycle multiplier over all simple cycles in 
the graph T. Then T is box expansive (hence f is hyperbolic) if and only i/M > 1. 

Thus if some F is box expansive, then each cycle multiplier for F is greater than 
one. However, some boxes in a cycle could have multiplier less than one, if others 
are large enough to compensate. In this case, in the euclidean metric, the map is not 
expansive along every edge. The handicaps are designed to spread the expansion 
out along the cycles, so that in the box metric, the map is expansive by at least L 
on every edge in the graph. This is why we use the term "handicap" . 

According to Corollarv ll.81 we can show whether / is box expansive simply by 
computing M. There do exist efficient algorithms for finding such a minimum, see 
15. However, in order to have explicit information about the hyperbolic structure, 
we still want to find a viable box expansion constant L < L, and build an expanded 
box metric, which we can do using the Handicap Hedging Algorithm. 

One weakness of this algorithm is that the box expansion constant, L > I, must 
be inputted in advance. Proposition 11.71 shows the ideal box expansion constant 
is the smallest average cycle multiplier, L. However, in Section Ffi. II we describe 
counterexamples which suggest that the only algorithm for explicitly computing 
£ is exponential (thus too inefficient for our examples). But then we describe an 
efficient method for finding a good approximation to the ideal L, in Section 1^^ 

Our approach for testing hyperbolicity has some similarities to work of Osipcnko 
( |17l I18| ). For / a diffeomorphism of a compact Ricmannian manifold M C K", 
he uses a symbolic image F of /, and develops a general algorithm for describing 
the expansion of / by approximating Lyapunov exponents and the Morse spectrum 
of the chain recurrent set. His algorithm can also be used to verify hyperbolicity. 
Our work here differs in that we are interested in developing and implementing 
efficient algorithms for families of polynomial maps of C (of degree greater than 
one). For this study, we found it more efficient to get expansion information and a 
hyperbolicity test by constructing a metric expanded by /. Osipenko's hyperbolicity 
test has the same computational complexity as finding the smallest average cycle 
multiplier for all simple cycles in the graph F (see Section IHTT) . 

In implementation, we control round-off error using interval arithmetic (lA). This 
method was recommended by Warwick Tucker, who used it in his recent computer 
proof that the Lorenz differential equation has the conjectured geometry (|22|)- In 
designing our algorithms, we must keep in mind the workings of lA. We thus give 
a brief description of lA in Section O 

To summarize the organization of the remaining sections: in Section |2] we show 
box expansion implies the standard definition of expansion, to prove Theorem 1 1.41 
in Section El we briefly describe interval arithmetic; in Section 01 we give our ba- 
sic algorithm, the Handicap Hedging Algorithm, for attempting to establish box 
expansion; in Section we obtain dynamical information from either success or 
failure of the Handicap Hedging Algorithm, proving Theorem ll.5l Proposition ll.71 
and Corollarv ll.81 in ScctionElwe compare and contrast ideal versus efficient meth- 
ods for determining a good expansion constant L; and in Section Owe discuss our 
implementation of the algorithm and give examples of output. 
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2. Box EXPANSION IMPLIES CONTINUOUS EXPANSION 

In this section, wc show box expansion iniphcs the standard definition of expan- 
sion. Throughout, let / denote a polynomial map of C of degree d > 1, and as in 
Definition ll.il let F be a box chain model of / on J, with vertex set V = {Bk}^^i, 
composing the box Julia set 23 = U^^^Bk D J. 

First, it is more natural for computer calculations, and reduces round-off error, 
to consider vectors in M^, rather than C, and use the L°° metric of M^, rather than 
euclidean. Hence, we consider 

\z\ =max{|Re(z)|,|Im(z)|}. 

Also, let ^^(5*, r) denote the open r-neighborhood about the set S in the metric in- 
duced by the above. This metric is uniformly equivalent to the euclidean metric |-|g, 
since -^ \-\^ < \-\ < \-\^. Thus neighborhoods are slightly different, but the topology 
generated by them is exactly the same, so they can nearly be used interchangeably. 
To prove Theorem 11.41 we use a partition of unity to smooth out a box metric 
<I>. First, we need a lemma from [H]. 

Lemma 2.1 ([HI)- There exists an rj > so that if Bk,Bj G V, with z E ?^(i?fc,7/) 
and f{z) g 3^(Bj, rj), then there is an edge from Bk to Bj in F. 

To prove this lemma, we used the assumption that / was a polynomial of degree 
d > 1, and the fact that by Definition ll.il there is a 5 > such that there is an edge 
from Bk to Bj if a ^-neighborhood of f{Bk) intersects Bj. Now, we can obtain: 

Lemma 2.2. Suppose F is box expansive. Then there exists a r > such that if 
Bk,Bj eV, and z e "Bn N(Sfc,T) with f{z) e J<{B.j,t), then for any v e T^C, 
\DJiv)\^>L\^\,. 

Proof. Among other requirements below, let r > be less than 77 from Lemma l2. II 
Then for z satisfying the hypotheses, there is an edge from Bk to Bj in F. 

Note since we are working in one dimension, Dxf = fix), hence box expansion 
yields that for x € Bk, (fij |/'(a;)| |v| > Lipk |v|, thus simply ipj |/'(a;)| > Lipk. 

Since CB is compact, V is finite, and f'{x) is continuous, there is an a > so 
that: 

(1) a = mm{ifij \f'{x)\ - Lifk- x e Bk, {k,j) e £}, 

(2) if T < 77 is sufficiently small, then for any j, \x — z\ < r implies that 
ip, If'ix) - nz)\ < a. 

Now z is not necessarily in Bk , but z e S , so suppose z G Bm and x e Bm H Bk 
such that |a; — z| < r. Then (pj |/'(a:;) — f'{z)\ < a; further, there is an edge {k,j) e 
£, hence x satisfies ipj \f'{x)\ — Lipk > «■ Combining these gives ipj \f'{z)\ > Lipk- 
Thus|Z)./(v)|. >L|v|,. D 

Now we use the r-overlap to convert a box metric into a riemannian metric. 

Definition 2.3. Suppose F is box expansive for a box metric $. Let r > be 
as given by Lemma 12.21 Define a partition of unity on S (F) by choosing smooth 
functions pk '■ C ^ [0, 1], for each box Bk G V, such that supp(pfc) C yi{Bk,T) and 
S/c Pkix) = 1, for any x € CB. Define the smooth function p = p($) : 53 ^ [0, 1] by 
p{x) = X]fe Pkix)(pk- Then p induces a riemannian metric on T23, with a smoothly 



6 S. L. HRUSKA 

varying norm, |-| , where if a; G 23 and w Cz TxTi then 

W\p = |v| pix) = |v| Y^ Pk{x)ipk = 51 ^fe(^) l^lfc • 

k k 

Note I -I and |-|^ are very close. They only differ in the small r- neighborhoods 
of the box boundaries, where the p metric smooths out the $ metric. 

It is also straightforward to show that both the p metric and the $ metric are 
uniformly equivalent to euclidean, with 

[mhii^k}^ M < Mp,|.U < (^maxj^fc}^ \-\. 

We establish Theorem 1 1 . 41 bv showing that if T is box expansive for a box metric 
$, then / is expansive in S for the metric p{^). 

Proof of Theorem \1.4\ Let |-L , t be as in Definition 12. 31 

Thus T is small enough that if a; e supp(/9fc), then for any j such that f{x) g 
supp(pj), we have ipj |Dj./(v)| > Lipk |v| , for any v G T^C Then if we set 

(fix = maxjyjfc : x S supp(pfc)}, and ipf^^ = min{.^_, : f{x) S supp(pj)}, 

we know ipf_x \Dxf(y)\ > Lipx |v| , for any v G T^C. 

Now if we use that ^j, pk{x) = J2j Pj{f{x)) = 1, we get the result: 

i 

> Lipx\v\>L\v\J2pk{x)ipk=L\v\p. 

k 

D 

Thus if any box model T of / on J is box expansive by some L > \, then / is 
expanding on J in the ricmannian metric p, and hence / is hyperbolic. 

3. Interval Arithmetic 

Interval arithmetic (lA) provides a natural and efficient method for manipulating 
boxes, and also for maintaining rigor in computations. The basic objects of lA are 
closed intervals, [a] = [a, a] G IK, with end points in some fixed field, K. An 
arithmetic operation on two intervals produces a resulting interval which contains 
the real answer. For example, 

[a] + [b] := [a + b,a + b\ 

[a] - [b] := [a - 5, a - b\ 

Multiplication and division can also be defined in lA. 

Computer arithmetic is performed not with real numbers, but rather in the finite 
space F of numbers representable by binary floating point numbers of a certain finite 
length. For example, since the number 0.1 is not a dyadic rational, it has an infinite 
binary expansion, so is not in F. 

Since an arithmetical operation on two numbers in F may not have a result in F, 
in order to implement rigorous lA we must round outward the result of any interval 
arithmetic operation, e.g. for [a], \b] £ IF, 

[a] + [b]:=[ia + bi,]a + b]], 
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where J,a;| denotes the largest number in F that is strictly less than x (i.e., x 
rounded down), and ]x] denotes the smallest number in F that is strictly greater 
than X (i.e., x rounded up). This is called lA with directed rounding. 

For any a; G M, let Hu11(j:) be the smallest interval in F which contains x. That is, 
if X G F, then Hull(a;) denotes [x,x]. If a; e R \ F, then Hull(a;) denotes [ja-j , ta-'T]- 

In higher dimensions, lA operations can be carried out component-wise, on in- 
terval vectors. Note a box in C = R^ is simply an interval vector. 

In designing our algorithms, every arithmetical calculation is carried out with 
lA, and thus we must think carefully about how to use lA in each situation. Our 
extensive use of boxes is designed to make lA calculations natural. However, lA 
tends to create problems with propagating increasingly large error bounds, if not 
handled carefully. For example, iterating a polynomial map on an interval which is 
close to the Julia set, J, can produce a tremendously large interval after only a few 
iterates (due to the expanding behavior of the map near J). Thus, in the remainder 
of the paper, after describing each algorithm we note how lA is being used. ^ The 
interested reader can find an abundance of materials on I A, for example [HI 1181 IT^ . 

4. The Handicap Hedging Algorithm 

In this section we describe our basic algorithm, the Handicap Hedging Algorithm, 
for attempting to build an expanded box metric for a box chain model F. 

The problem of finding a set of handicaps {(fik} defining a box metric ^ for which 
F is box expansive by a given L > 1 is strictly a graph theoretic problem. We want 
to find handicaps so that for every edge {k,j) G £, for every point z G Bk, and for 
every vector v G T^Bk, we have 

^,|i?,/(v)H(^,|D,/||v|>L^fc|v|, 

or equivalently, if A^ = min^gSj, |/'(^)| is the multiplier of Bk, we want 

(1) ipj > L(pk/>^k- 

For polynomial maps of C, a lower bound for the multiplier Afc in a box is easily 
calculated with interval arithmetic. After the Afe are calculated, we can forget the 
map / and simply think of F as a strongly connected directed graph endowed with 
edge weights, S,k,j = L/Xk, if {k,j) € £. In the case of maps in one dimension, 
the edge weights are all equal along all edges emanating from a single vertex, i.e., 
Cfc.i = ?fc,2 = • ■ ■Ck,n- However we describe the algorithm for the more general 
situation, where ^k.j j^ ^k,m, since it is not more difficult, and it is useful for higher 
dimensional applications (for example, in jS]). 

Definition 4.1. Fet F be a directed graph with positive edge weights S ~ {£,k,j '■ 
{k,j) e £}. Fet $ = {(fik- Vk e V} be a set of positive vertex weights, called 
handicaps. We call $ consistent handicaps (for F, S) if ipj > ^k j fk, for every edge 
{k,j)G£. 

If equality holds on every edge, then we call $ strict handicaps. 

Then F is box expansive if there exists a set of consistent handicaps <& for F, 
given the edge weights suggested by ^. 

To attempt to find consistent handicaps given F and the edge weights S, we 
break the problem into a finite induction. Step consists of finding an initial set 



All of our lA computations use the PROFIL/BIAS package, available at I2UI . 
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of strict handicaps for some spanning tree, To, of T. At steps n > 0, we choose a 
graph r„ such that To C r„ C F, then seek a consistent set of handicaps for r„. 

Definition 4.2. A directed graph Fq is an arborescence if there is a root vertex vq 
so that for any other vertex v, there is a unique simple path from vq to v. Such a 
graph is a tree, and must have exactly one incoming edge for each vertex v ^ vq. 

If F is strongly connected, then for each vertex vq in F, there is a minimum span- 
ning tree Fq with root vertex vq which is an arborescence (simply perform a depth 
first or breadth first search from vq)- We call such a Fq a spanning arborescence. 

Definition 4.3. Let F be a finite, strongly connected directed graph. Let Fq C 
Fi C • • • C Fjv = F be a nested sequence of subgraphs of F such that Fq is a 
spanning arborescence of F, and for n > 1, V(F„) = V(F), and £(F„) is formed by 
adding one edge of £(F) \ £(F„_i) to £(F„_i). We call such a sequence an edge 
exhaustion of F. 

Note by the above definitions that an edge exhaustion exists for any finite, 
strongly connected directed graph F. Note also that since Fq is a spanning ar- 
borescence, each F„ is edge connected. 

Algorithm 4.4 (Recursively hedging handicaps via an edge exhaustion). Let Fq C 
Fi C • • • C Ftv = F be an edge exhaustion of F. 

Base Case: Inductively construct a set of strict handicaps, $" = {f^} for Fq, 
by choosing any i^fj for the root vertex vq, then pushing this value across Fq by 
multiplication with successive edge weights. That is, if w^ e V(Fo), and v ^ vq, 
then since Fq is an arborescence, the predecessor, ^^(fc), is uniquely defined such 
that the edge {■^{k), k) is in Fq. Then set vj° = ^k,-K{k) V^ik)- 

For our box expansion application, we pick a random vertex as the root, and 
start with ip^ ~ 1, i.e., take the euclidean metric on that box. 

Inductive Step: Suppose we have a set of consistent handicaps (f)""^ = {ip^~~^} 
for F„_i. We attempt to define a set of consistent handicaps <&" = {^p^} for F„, by 
adjusting the set $"~^, using a process we call recursively hedging the handicaps. 
First, start with the set of temporary handicaps <&" := $"^1 on F„. Let {va,Vb) 
denote the edge in F„ that is not in F„_i (this edge is unique by definition of edge 
exhaustion). Then we adjust the 6th handicap, by doing 

set V^b = la,b (/j" 
If the second line above is performed, i.e. p^ is strictly increased, then we say 
we have hedged ip^ along the edge {va,Vb). But now if (p^ is hedged, then in order 
to keep consistency of the handicaps, for each vertex Vc in F„ which is adjacent to 
Wb, we may need to hedge (p" along the edge {vb,Vc)'. 

set p^ = Ub ^a 

for each Vc in F„ and in the adjacency list of ff, do 
if Vc < ib,c H=l 

set p^^ = ib,c Vb 

The hedging process becomes recursive here, for if any (/?" is increased, then we 
may need to hedge handicaps along each edge in F„ emanating from Vc, etc. The 
worst case scenario is that every vertex reachable from Vb along every path in F„ 
would need to be hedged. 
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However if at any step in the hedging, we find a vertex Vx whose handicap need 
not be increased, i.e., ip^ > £,'k{x),x Vwrr) aheady, then we may stop, since the 
handicaps of vertices reachable from v.^ will not need to be increased. This greatly 
saves on computational time. Thus, our process so far is: 
Hedge(a, b) = 

if Vb < Ca,6 ^a 

set V3^ = Ub (fa 

for each Vc in r„ and in the adjacency list of ff, do 
Hedge(6, c) 

Will this recursive procedure always terminate in a consistent set of handicaps 
<I>" for r„? Of course not, since it is possible that there does not exist any set of 
consistent handicaps for T with edge weights S. There is one possible obstruction: 
suppose that for the new edge in r„, edge {va, Vb), we travel away from Vb, hedging 
the handicaps along every edge in some path, Va,Vb,Vc, . ■ . ,Vx, and then discover Va 
in the adjacency list of v^. This means every handicap on a path in r„ from Vb to 
Vx has been increased, since we stopped the search if it was not. The handicap at 
Va cannot be increased, since then we would have to hedge ip'^ along (va,Vb) again, 
and an infinite loop of hedging along this path through Vx would occur. 

Thus the cycle {va, Vb, Vc, ■ . ■ , Vx,Va} could be an obstruction to consistent hand- 
icaps, for upon seeing Va from Vx, we can only check if (/?" > ^x,a 'fix already. This 
leads us to the following algorithm for the n*'* inductive step: 

do Hedge(a, b, a), where 
Hedge(u, w, a) = 

if ^v > S,u,v ^ui f^^cn return 1 
else if (w = a), then return 
else 

set V3^ = ^u^y (^2 

for each lu in r„ and in adjacency list of v do 
Hedge(w, w,a) 

Thus the n* step is terminated when either a cyclic obstruction is found (at 
the "return 0" line above), or when all necessary hedgings have successfully been 
performed (if Hedge(a, 6,a) returns 1). If no cyclic obstructions are found, then 
$" is a consistent set of handicaps. 

Thus Algorithm 14. 41 either terminates in a cyclic obstruction at some r„, or pro- 
duces consistent handicaps for all of F = Fw . This dichotomy leads to Theorem ll.5l 
proved in Section [S] 

The Handicap Hedging Algorithm consists simply of applying Algorithm 14.41 to 
the case of showing box expansion for a box chain model F of / on J. We describe 
this process more explicitly with the following pseduo-code. Comments are paren- 
thetical and to the right. The routine consists of a main function, BuildMetric, 
and its two helper sub-routines, SpanTree and Hedge. SpanTree creates an ar- 
borescence using a breadth first search and a queue (first-in = first-out) to traverse 
the graph. Hedge uses a depth first style to check the rest of the edges. See 0] for 
background on graph searches, and note the remark following the pseudo-code on 
why these styles were chosen. 

Algorithm 4.5 (The Handicap Hedging Algorithm). 
BuildMetric(r,L): 
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for every vertex u in V 
set color [it] to white 

for each vertex v in adjacency hst of u do 
set check [w][v] = 
set 1^0 = 1 

do SpanTree(vcrtex of F) 
for every edge (a, h) in T 

if (check[a][6]= 0) then [edge {a,b) not in previous step of exhaustion) 

set check[a][6] = 1 (put (a, 6) into current step of exhaustion) 

if (Hedge(a, b, a, L) = 0) then return 
return 1 {^f get this far, then successful) 

SpanTree(M): 

put vertex u on the queue, Q (at the end) 
set color [u] to gray 
while Q 7^ 

let V be the head (first element) of the queue 
for each vertex w in adjacency list of v do 
if color [to] is white 

put vertex w on the queue, Q (at the end) 
set color [w] to gray 

set check[t)][w] = 1 (record edge {v,w) is in spanning tree) 

set ifw = Ltp^/Xy 
remove first element from queue (remove v) 
set color [w] to black 

Hedge(u, w, a, L): 

if (y^v > Lifu/Xu) then return 1 (i.e., edge (u,v) is already ok) 

else if (v = a) then return (we cannot increase (fa, so fails) 

else 

set Lpv = Lip.u/\u (increase ipv) 

for each w in adjacency list of v do 

if (check[v][w] = 1) then ((v,w) is in current step of exhaustion) 
if (Hedge(w, w, a, L) — 0) then return 
return 1 (if get this far, then successful) 



Remark 2. In Span Tree we maximize efficiency by using a breadth first style 
search to traverse the graph, instead of depth first. Running depth first search on 
a typical box chain model T tends to produces a spanning arborescence with a very 
long path, whereas breadth first search constructs an arborescence with paths of 
minimum length. This is because for a polynomial map /, the typical T has a large 
number of vertices and a small bound on the out-degree of vertices (related to the 
bound on the derivative of the map in 'B). In addition, F is strongly connected. 

Long paths creates two problems. First, creating such a tree can causes a memory 
overflow, due to many nested recursive function calls. Additionally, even when there 
is no crash, a large path tends to create an initial metric farther from Euclidean. 
See Example 16. II in Section IHl for an exploration of the latter phenomena. 
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Recall in our implementation, wc control rounding using interval arithmetic (Sec- 
tion I^J. In the above, we round up to ensure the inequality is satisfied, i.e., check 
ifiy > sup(Hull(L) * Hull((p„)/Hun(A„)). 

5. Characterization of box expansion 

Recall from Dcfinition ll .6l that the multiplier of a box Bk is Xk = min,gBj. | f'{z) \ , 
and for an n-cycle of boxes {Bq -^ . . . ^ B„-i -^ Bn = So) hi F, we call 
(Ao • • • A„_i) the cycle multiplier and (Aq • • • A„_i)^/" the average cycle multiplier. 
Also, a cycle in a graph is called simple if is it composed of distinct vertices. 

We can now specify the implications of success or failure in the Handicap Hedging 
Algorithm, proving Theorem 1 1.51 ProDOsition ll.71 and Corollarv ll.8l 



Proposition 5.1. For L > 0, the Handicap Hedging Algorithm, i.e., Algorithm \4.5\ 
or Algorithm \4.4\ "with edge weights S = {^k,j = L/Xk : {k,.j) G £}, either constructs 
a set of consistent handicaps showing T is box expansive by L, or finds an n-cycle 
of boxes with cycle multiplier less than L" . 

Proof. We observed in the description of Algorithm 14.41 that the only obstruction 
to building consistent handicaps is if in adding an edge (u,w), we find a cycle of 
boxes {u = Bq ^ V = Bi ^ . . . ^ Bn-i -^ Bn = Bq), such that, holding (pQ fixed, 
in order to keep consistency the metric handicaps must be increased along every 
edge in the cycle. That is, (pk+i = Lipk/Xk for < k < n — 2, and we have the 
failure (pQ < L(^„_i/A„_i. But then 

Lip„-l L'^'Pn-2 L"ipQ 

•fo < —. 



-1 A„_iAn-2 An-lA„_2 ■ • ■ Ao 

Hence, A„_iA„_2 • • • Ao < L". D 

Lemma 5.2. Let {Bq —>■...—>■ Bn-i -^ Bn = Bq) be an n-cycle of boxes in T , 
such that consistent handicaps {ipQ, . . . , ipn-i, fn = ^o} can be chosen to show f 
box expands by L > I along the cycle. Then its cycle multiplier is at least L". 

Proof. By hypothesis, wc know ipk+iXk > Lipk, k € {0, . . . , n — 1}. Thus, 
-^0 ■ ■ ■ -^)i-i ^ L — L — ■ ■ ■ L- 



^1 ^2 P>n = VQ 

Cross cancellation and simplifying leaves only Ao • • • A„_i > L". D 

Lemma 5.3. For L > 0, there is a cycle of boxes in T with average multiplier less 
than L if and only if there is no box metric which T expands by L. 

Proof. First, Lemma l5 . 21 gives the forward implication immediately; for, if there is 
a box metric which F expands by L, then every cycle must have average multiplier 
at least L. Conversely, suppose there is no box metric which F expands by L. Then 
by Proposition 15. II there exists a cycle with average multiplier less than L. D 

Proof of Theorem \1.5l Proposition l5 . ll shows that given L > 0, the Handicap Hedg- 
ing Algorithm either constructs a box metric for which F is expansive by L, or finds 
a cycle with average multiplier less than L. Lemma l5.3l shows the latter is equivalent 
to saying there is no box metric which F expands by L. D 
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Proof of Proposition \1.7\ This follows directly from Lemma 15.31 Indeed, let L > 
be given. First, assume L < L. Then all cycles have average multiplier at least L, 
so by Lemma |5. 31 there exists a box metric which F expands by L. Next, assume 
L > L. Then the cycle of boxes with average multiplier Z has average multiplier 
less than i, hence by Lemma 15.31 there is no box metric which F expands by L. D 

Proof of Corollarv \l.y[ This follows easily from Proposition ll.71 using the fact that 
a > 1 if and only if a^/™ > 1, for any positive a and m. Let H denote the smallest 
average cycle multiplier over all simple cycles in the graph. By Proposition 11.71 
F is box expansive if and only if £ > 1. Let m be the length of the cycle with 
multiplier M. Let n be the length of the cycle with average multiplier £. Then 
M^/™ > £ and M < L". Suppose F is box expansive. Then M^/" > ii > I, hence 
M > 1. Conversely, suppose M > 1. Then £," > M > 1, hence £, > 1. Thus F is 
box expansive. D 

6. Finding a good box expansion constant L 

In theory and in practice, a box expansion constant more appropriate for a 
particular collection of boxes yields a "better" metric. In particular, in running the 
program Hypatia it is easy to find behavior like the following. Suppose for some 
map, we successfully build a metric with box expansion constant L = 1.2, but the 
resulting handicaps range between 4.3 x 10~^* and 0.6. Depending on the map /, 
trying L = 1.5 could improve the handicaps to be in the range 0.068 to 0.55. This 
yields a metric more computationally tractable, and closer to euclidean. Why do 
we see this behavior? Consider the example: 

Example 6.1. Suppose there is a path of boxes Bq ^ . . . ^ _B„ such that (for 
simplicity) all multipliers are approximately the same value A. Then suppose we 
try to show box expansion by some L < X. We would define the handicaps in the 
boxes of this cycle so that (pk+i > Lipk/X, for < fc < n — 1. Thus we would get 



fi = — ^, f2 = —r- = — r^-7 ■■■,'Pn 



A '"^ A _^2 '■■■'-" x^' 

But since L < X, large n leads to L" <C A" and (/?„ <^ ipo. Thus, if we use an L 
which is "too small", then the handicaps plummet. 

In fact in running Hypatia, we have found that it is easy for a bad choice of L to 
lead to handicaps which are so large or small that the machine cannot distinguish 
them from or oo. Checks must be put in place in the algorithms to flag such 
occurrences. This can be guarded against somewhat by using breadth first search 
style algorithms as much as possible, rather than depth first search. In particular, in 
SpanTree, if a depth first search is used, then the spanning arborescenccs usually 
contain very long paths. On longer paths, as Example l6.1l illustrates. the handicaps 
are more likely to get unmanageable. 

6.1. An optimal, yet impractical solution. Proposition 11.71 implies that the 
ideal box expansion constant, H, is equal to the smallest average cycle multiplier 
over all simple cycles. This can theoretically be computed, since the graph is finite. 
However, in this section we shall see the computation seems unwieldy. 

For a first approach to finding iL, recall that there do exist efficient algorithms 
for finding the cycle with the smallest cycle multiplier (0]). Unfortunately, the 
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Figure 1. Graph illustrating Examples 16. 21 and 16.41 

following example shows that the smallest average cycle multiplier cannot be simply 
derived from examining the cycle with the smallest cycle multiplier. 

Example 6.2. Consider the graph with three vertices {0, 1, 2}, and with the four 
edges {(0, 1), (1,2), (0,2), (2,0)}, shown in Figure [2 It has two cycles: 0^1^ 
2 ^ and -> 2 ^ 0, with multiphcrs: Aq = 1/2, Ai = 3/2, A2 == 8. The cycle 
multipliers are: (A0A2) = 4 and (A0A1A2) = 6. Thus the smallest cycle multiplier 
is along the cycle 0-^2^0. The average along this cycle is 4^/^ = 2. However, 
the average cycle multiplier along the other cycle is (AoAiA2)^/'^ = 6^/'^ = 1.82. 

Thus, we need an algorithm for computing the smallest average cycle multiplier 
from scratch. For our second approach to this computation, we might posit that we 
could adapt the algorithm used to compute the smallest cycle multiplier, to instead 
compute the smallest average cycle multiplier. This algorithm involves considering 
"best paths" during a breadth first search of the graph. However, in the following 
we explain how the "average" prevents this shortcut from working. 

Definition 6.3. If P = {Bq -^ . . . ^ Bn) is a path in the graph F, then the path 
multiplier IS the product of the multipliers along the path, \{P) := Aq ■ ■ • A„_i, and 
the average path multiplier is (A(P))^'". 

Let Bi^ and Bj be vertices in a graph F. Let P^. j be the path from Pj. to Bj 
with the smallest average path multiplier of all paths from Bk to Bj . The following 
example shows that there exist graphs F such that the cycle containing Bk and Bj 
with the smallest average multiplier docs not contain the path Pk,j- 

Example 6.4. Again consider the three vertex graph in Figure Q] used in the 
proof of Example 16.21 It has two cycles: ^ 1 ^ 2 ^ and ^ 2 ^ 0, and has 
multipliers: Aq = 1/2, Ai = 3/2, A2 = 8. Comparing paths from vertex to vertex 
2, we see Ao = 1/2 and (AqAi)^/^ = V3/2. The smallest average path multiplier 
is along the path 0^2. However, Example 16.21 shows the smallest average cycle 
multiplier is along the cycle ^ 1 ^ 2 ^ 0. 

Example 16.41 suggests that an algorithm to compute the smallest average cycle 
multiplier in a graph must compute the average cycle multiplier for each simple 
cycle separately. A simple combinatorial argument shows that the number of simple 
cycles in a graph (even a sparse graph) is exponential in the number of vertices of 
the graph. The graphs created for polynomial maps of C are large enough to make 
an exponential algorithm prohibitively inefficient (sec Section|Z|for several examples 
of such graphs). Thus, we cannot expect an efficient algorithm for determining the 
ideal box expansion constant. 
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The previous exploration can also be used to examine the computational com- 
plexity of Osipenko's method Ql7lll8) 'l for approximating Lyapunov exponents, or 
the Morse spectrum, on J. and the subsequent hyperbolicity test. If C = i?o — > 
. . . ^ B„-i — > Bn ~ Bq is an n-cycle of boxes in F, then the characteristic Lya- 
punov exponent of the cycle is bounded below by 

1 "~^ 1 f"^^ \ 1 

— 7 In(Afc) = — In I I A^ = — In (cycle multiplier of C) . 

n ^-^ n \ -'■-'- / n 

fe=o \fe=o / 

Approximating the Morse spectrum of / means computing the characteristic Lya- 
punov exponent of each simple cycle in F, and since there are an exponential number 
of simple cycles, this is an exponential-time algorithm. Osipenko's test for hyper- 
bolicity says that if the minimum exponent over all simple cycles is positive, then 
the map is hyperbolic. However, the graph of Example 16 . 41 also suggests there is no 
shortcut to efficiently computing the minimum exponent, since taking the log and 
dividing by the cycle length for this example produces the exact same complications 
as taking the n*^ root (the specific numbers can be easily checked). 

6.2. Approximate, efficient solutions. In order to efficiently determine a suc- 
cessful box expansion constant L, we need a good starting guess. We have a lower 
bound of 1 for L. CoroUarv 11.81 gives a way to get an upper bound: if M > 1 is 
the smallest cycle multiplier, realized by a cycle of length m, then the ideal box 
expansion constant (i.e., the smallest average multiplier) £ satisfies 1 < L < M^/™. 
Alternatively, we can get an upper bound for L in a very simple way. For a 
polynomial map /, the Lyapunov exponent A measures the rate of growth of tangent 
vectors to J under iteration. A description of the one- variable case is given in |21| . 

Theorem 6.5 ([21 E3)- ^o*" "^ polynomial map f : C ^ C of degree d > 1, the 
Lyapunov exponent A satisfies A > logd, with equality if and only if J is connected. 

Thus for degree d maps with connected J, L = d is an upper bound. 

One straightforward way to obtain a good L in a preset number of steps is basic 
bisection. Keep track of loL, the most recent working L, and hiL, the most recent 
failing L. Lower L halfway to loL when box expansion fails, and raise it halfway 
to hii when it succeeds. Start with \oL = 1 and hiL = d or M^/™ if computed. 

An alternative to straight bisection is to utilize the information that Algo- 
rithm 14.51 BuildMetric, already discovers in a test for box expansion. Indeed, 
if box expansion fails for some L, then we realized in Proposition II . 71 that it is due 
to a cycle, hadcycle, with average multiplier L' less than L. But we can easily adapt 
the algorithm to compute and return the average multiplier of hadcycle, L' . Then 
if i' < 1, we know the map is not box expansive on F and we can stop. Otherwise, 
on the next pass instead of lowering by some arbitrary amount we may simply test 
by this new average multiplier L' . Better yet, we can test by the minimum of L' 
and L minus some preset step size (5, in order to prevent increasingly small steps 
down. We refer to the above procedure as the algorithm Find-L-Cycles. 

Proposition 6.6. Let L be the minimum average multiplier over all simple cycles 
in the graph F. // 2 > iL > 1 + (5, then Find-L-Cycles shows box expansion by 
some L within 6 of L, in at most 1/5 trials of BuildMetric. 

Proof. Since we are looking for an expansion amount in the interval [1,2], and 
decrease by at least 5 at each step, we perform at most 1/5 attempts. Suppose 
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one of these attempts is successful. That is, suppose F fails to box expand by 
some Lo < 2, and thus outputs an average cycle multiplier of Li < Lq. Since Z 
is the minimum, H < Li. Suppose Hypatia verifies successfully box expansion by 
L = minjLi, Lq — S}. Then by Lemma 15.21 we have L > L. Thus we have: 

Lq > Li > £j > L ~ minjLi, Lq — S}. 

Thus, either we were lucky and L = Li and we have shown box expansion by 
exactly that amount, ox Lq > Li > L > Lq — 5 and we have shown expansion by 
Li within 5 oi L. D 

The best way to control round-off error in the above with interval arithmetic (see 
Sectional) is to "round down" , since a box expansion constant larger than the aver- 
age cycle multiplier would fail. That is, set Li = Inf((Hull(Ao) • • • Hull(A„_i))^/"). 

7. Examples of running Hypatia for polynomial maps 

In this section, wc describe the results of applying the algorithms of this paper to 
some examples for quadratic and cubic polynomials, Pc{z) = z'^ + c and Pc,a{z) = 
z^ — 3a^z + c, using our implementation in the computer program Hypatia. In 
particular, in all of the examples of this section, we used the Handicap Hedging 
Algorithm and the method Find-L-Cycles f Section 16.2(1 . to try to produce an 
L > 1 and a box metric for which the constructed F was box expansive. Table^^t 
the end of the section summarizes the data for all the examples. 

All of the computations described in this section were run on a Sun Enterprise 
E3500 server with 4 processors, each 400MHz UltraSPARC (though the multipro- 
cessor was not used) and 4 GB of RAM. ^ 

7.1. Producing the box model F. First we summarize how we used the box 
chain construction from [H| to produce box models F for these maps. The box 
chain construction is in fact an iterative process. We began by defining a large box 
23o in C such that J C Sq. Then for some n > 1, we place a 2" x 2" grid of boxes on 
"Bq. The construction then builds a graph F consisting of a subcoUcction of these 
grid boxes, which is a box model of J, according to Definition ll.il If the boxes are 
sufficiently small, then every invariant set disjoint from J is disjoint from !B(r) (for 
example, attracting periodic orbits). Further, we can produce an improvement of 
F by subdividing all (or some) of the boxes in F, then repeating the construction 
to produce a new graph F' such that S(F') C 'B(F). 

For quadratic polynomials, it is easy to check that if \c\ < 2, then the filled Julia 
set is strictly contained in the box [—2,2]^. For the cubic polynomials Pc^a{z) = 
z^ — ia^z + c, one can also check that the filled Julia set is contained in [—2,2]^ 
whenever \c\ < 2 and \a\ < (2/3)^/^. We denote ©2 = [-2,2]^. In each of our 
examples, we began the box chain construction with some 2" x 2" grid on CBq = ID)2- 

7.2. Selective subdivision. We also introduce here two modifications to the basic 
box chain construction, in order to improve efficiency. The key idea for both is 
selective subdivision, i.e., rather than subdividing all of the boxes to form the next 
level, we only subdivide the boxes where the dynamics is behaving badly. 

Thus we must determine which boxes are most obstructing the construction of 
a hyperbolic metric. One approach is to choose boxes which are closest to a sink 



The server was obtained by the Cornell University mathematics department through an NSF 
SCREMS grant. 
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Figure 2. A box Julia set for S(r), for the map z ^^ z^ - 1.755, 
and boxes from a 2^-"^ x 2^^ grid on ©2- T is box expansive. 

cycle, say boxes which have bounded midpoints after several images, and subdivide 
only those boxes. We call this procedure sink basin subdivision. See Example 17.41 
Alternatively, the Find-L-Cycles algorithm, in addition to improving the metric 
on the given graph, can also be used to determine a selective subdivision procedure, 
which we call weak cycle subdivision. After Find-L-Cycles is run, we can easily 
identify the boxes involved in the cycle with the weakest cycle multiplier. Then 
subdividing these boxes should help improve the expansion amounts. We utilize 
weak cycle subdivision in Examples 17. 21 and 17. 31 

7.3. Examples for quadratic polynomials. The quadratic polynomial Pc{z) = 
z^ + c has one critical point, at zero, thus has at most one (finite) attracting cycle. 

Example 7.1. The aeroplane, Pc for c = —1.755, has a period 3 sink. Shown in 
Figure[2|is a box Julia set 'B(r) consisting of boxes from 2^^ x 2^^ grid on D2. This 
box chain model V has 43,000 vertices and 260,000 edges. T is box expansive, with 
box expansion constant L = 1.05069. The computation took less than 200 MB of 
RAM and 6.5 minutes. 

Example 7.2. The quadratic polynomial P{z) = z-^ — 1 is called the basilica. This 
map has a period 2 attracting cycle <-^ — 1. The picture on the upper left of 
Figure 13 is a heuristic sketch of J, drawn using the program Fractalasm (available 
at P). Shading is according to rate of escape to infinity, or to the two-cycle. 

The picture on the lower left of Figure 13 is a box Julia set S(r), composed of 
selected boxes from a 2'^ x 2'^ grid on ID2. F has 1,800 vertices and 14,000 edges. 
We could have easily made a finer picture, but this rough box Julia set is enough 
to prove hyperbolicity. Indeed, this box chain model F is box expansive, with box 
expansion constant L = 1.14067. The associated box metric has handicaps, ^pk^ 
with minimum 0.019, average 0.049, and maximum 1. This initial computation 
took only 8 MB of RAM and less than a minute of CPU time. 

For a deeper understanding of this metric, we created a picture of the box Julia 
set, with shading of each box according to the handicap for that box, shown on the 
lower left of Figure 13 Boxes with smaller handicaps are shaded darker. 

Example 7.3. A Cantor Julia set near the cusp of the Mandelbrot set is c = .35. 
The upper right of Figure |31 is a Fractalasm sketch of this Julia set. This is an 
example in which we can use weak cycle subdivision to show hyperbolicity more 
quickly. The box chain model from the 2^ x 2^ grid on D2 fails to be box expansive 
by L = 1 due to a bad cycle of length only 1. We had the program subdivide just 
that box, and again the resulting box chain model failed for L = 1 due to a length 
1 bad cycle. But, upon subdividing that one box, we found that the resulting box 
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Figure 3. Left: sketches of J for the map z i-^ z^ — 1; Right: 
sketches of the cantor J for the map z i— *■ z^ + .35; Upper: heuris- 
tic sketches of J, drawn with the program Fractalasm; Lower: a 
box-expansive box Juha set 23 (F) for each map, where shading is 
determined by handicaps, with smaUer handicaps colored darker. 
Note the asymmetry in metric shades. This is not dynamically 
significant, but an artifact of the tree construction for assigning 
handicaps. Boxes on the lower left have side length 4/2^, and 
boxes on the lower right have side length A/ 2^. 



chain model is box expansive (by L = 1.00778). On the lower right of Figure |21 is 
the hyperbolic box Julia set, with shading according to the value of the handicap 
in each box. Boxes with darker shading have handicaps closer to zero. 



7.4. Examples for cubic polynomials. The cubics Pa^c{z) 



3a^z + c have 



two critical points, at ±a, thus have at most two attracting cycles. When a = 0, 
there is a clear correspondence between z^ + c and z^ + c. 

Example 7.4. A seemingly rabbit-like cubic is Pa.c, c = —.44— .525i, a ~ .3i. This 
map has an attracting period three cycle. Shown in Figure 01 on the left is the 
box Julia set 'B(ra,c) with boxes from a 2^° x 2^^ grid on D2. However, this F 
is not hyperbolic because B contains one of the critical points, a ~ .3i. We then 



18 



S. L. HRUSKA 




Figure 4. Each image above is a box model of J and the period 
3 sink for a cubic polynomial Pa,c- Points heuristically found to 
be in the filled Julia set are shaded lighter, to illustrate J. Left: 
the box Julia set :B for c = -.44 - .525i,a = .3i, for a 2^" x 2^" 
grid on ©2 contains a critical point, so the box chain model F is 
not box expansive. Center: a refinement of the left picture, with 
boxes of side length 4/2^° and 4/2^^, is box expansive. Right: the 
box Julia set S for the map c = —.38125 + .40625*, a = 0.5i, with 
boxes of side length 4/2^ and 4/2^° is box expansive. 



performed sink basin subdivision and produced a box chain model F' which is box 
expansive. The latter box Julia set S' is shown in the center of Figure 0| 

Example 7.5. The cubic polynomial Pa,c, with c = — .38125 + . 40625z, a = Q.5i also 
has an attracting cycle of period three, but here the Julia set is disconnected. The 
most efficient method to get a hyperbolic box chain model F is to first subdivide 
all boxes uniformly to obtain a 2^ x 2^ grid on D2, and then perform sink basin 
subdivision. Shown on the right in Figure 0| is the resulting box Julia set. 
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Table 1. Data for the box chain models described in Section [7| Here the box depth n is the number such that the 
boxes are of size 4/2" (from a (2" x 2") grid on ©2)- 
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